Observing Nucleon Decay in Lead Perchlorate 
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Lead perchlorate, part of the OMNIS supernova neutrino detector, contains two nuclei, 208 Pb and 
35 CI, that might be used to study nucleon decay. Both would produce signatures that will make 
them especially useful for studying less- well-studied neutron decay modes, e.g., those in which only 
neutrinos are emitted. 

PACS numbers: 13.30.Ce, 14.20.Dh 



I. INTRODUCTION 

Studies of nucleon decay (see, e.g., 0, 0> 0, 0, E|) 
are among the most important in physics, in that they 
provide direct tests of fundamental theories of particle 
physics. The different particle theories make rather dif- 
ferent predictions as to what the nucleon decay half-lives 
might be, and even which decay modes would be expected 
to dominate. Although a number of searches for nucleon 
decay have been performed in the large detectors that 
exist, no convincing evidence has yet been presented for 
its existence. The best limits from these studies have in- 
volved decay modes in which protons decay to relativistic 
charged leptons (see, e.g., 0), either as individual pro- 
tons or as protons in nuclei, as the resulting Cherenkov 
radiation would produce definitive signatures, typically 
of order 10 33 y. However, nucleon decays in composite 
nuclei 0, 0, EJ might produce clear evidence for the 
existence of such effects that could not be obtained from 
decay of isolated protons. Furthermore, these might al- 
low extension of the existing limits of the decay branches 
in some instances, even with relatively small detectors. 

In this study we focus on the n — * v + v + v decay, 
which has a less well-known half-life because its decay 
can often be masked by backgrounds, and would in any 
event be difficult to detect via the means used to observe 
nucleon decay into relativistic leptons. It is not neces- 
sarily expected to be the dominant mode of decay, but 
it is the least-well-determined mode, so is the primary 
limitation to the mode-independent half-life. Moreover, 
there are non-standard models in which unusual decay 
modes may turn out to be the dominant ones (see, e.g., 
Illf). In some models the three neutrino decay mode can 
emerge naturally as the dominant one |12| . 

Decays of the three neutrino mode were searched for 
a decade ago, |13| . then their limit was improved to its 
current experimental value of 5xl0 26 y 0]. Recent sug- 
gestions for studying this decay mode have included sig- 
natures that would result from decay of a neutron in the 
O in H2O |9j of Super-Kamiokande and in the C in the 
CH2 [l3 of KamLAND. These suggestions both involve 
signals that would be generated by rather weak branches 



resulting from nucleon decay. 



II. SIGNALS FROM LEAD PERCHLORATE 

In this paper we study two nuclides, 208 Pb and 35 C1 
which, we show, would produce special responses to neu- 
tron decay. We present the signatures that might re- 
sult in the n — > v v v decay. We find that 208 Pb would 
have good sensitivity to this mode, but 35 C1, while not 
producing as strong a limit, would exhibit an unusually 
definitive signature. Both nuclei are part of a planned 
detector, lead perchlorate, LPC, Pb[C104]2, that will be 
part of OMNIS, the Observatory for Multiflavor Neutri- 
nos from Supernovae [l4|- LPC is a colorless liquid that is 
highly soluble in water. The properties of such a mixed 
liquid have been studied extensively [Tsj: it was found 
that the LPC will produce Cherenkov light from rela- 
tivistic particles, e.g., electrons and muons, that might 
be produced in many nucleon decays. In addition, each 
neutron produced by neutron emission from a nucleus 
within the detector will be captured on the 35 CI within 
tens of /is, producing 8.6 MeV of 7-rays. An energetic 
7-ray will also produce Cherenkov light, with a pattern 
that is indistinguishable from that from the relativistic 
leptons in LPC. 

In general we can write the following expression for the 
neutron decay lifetime, r„, as 



T n /Br > N n e n R det eie2, 



(1) 



where N„ is the number of neutrons in, e.g., 208 Pb or 
35 C1, e n is the fraction of those that can decay into de- 
tectable signatures, Rdet is the observed event rate, and 
ei and £2 are the detection efficiencies for the two (or 
more) signatures of the decay. The factor £3 accounts for 
the possibility of multiple neutron emission from decay 
of a neutron in 208 Pb, each with detection efficiency £2- 
The factor Br is the branching ratio for the decays that 
go to the specific decay mode being studied. In 208 Pb 
decay, this could, e.g., refer to the branching ratio for 
producing one neutron and a subsequent 7-ray with suf- 
ficient energy to be observed; we have assumed E ex < 
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3 MeV in the daughter nucleus for this criterion to be 
satisfied, as virtually every level with excitation en ergy 
above that will produce at least one 3 MeV 7-ray (16j |. 
In 35 CI decay, Br refers to a decay to states that will 
produce a 7-ray with at least 3 MeV of energy together 
with decay to the 34 C1 ground state. ( 34 C1 has an iso- 
meric state at 0.146 MeV, to which roughly half of the 
highly excited states will ultimately decay. However, it 
has a much longer half- life than the 1.5 s half-life of the 
ground state, so will be assumed not to be useful for the 
present discussion.) 

In order to estimate the probability of observing de- 
cays from 208 Pb and 35 C1 we will have to consider the 
probability of population of sufficiently highly excited 
states in the daughter nuclei: 207 Pb and 34 C1. Gen- 
erally the more deeply bound the nucleon that decays 
the larger will be the excitation energy in the daughter 
nuclide. The nuclide that decays will be a nuclide in 
the parent nucleus, so the resulting "state" of the ini- 
tial nucleus, 208 Pb or 35 CI, minus one neutron will map 
onto actual states in the daughter nuclide (assuming the 
neutron decay leaves the resulting nucleus intact; a rea- 
sonable assumption for the decay mode that produces 
three particles that interact only through the weak in- 
teraction). However, consider the energy E 'decay of the 
neutron decay products (e.g. the 3 neutrinos) when the 
final nucleus A-l is excited to the excitation energy E exc . 
Then E decay = M n — (S n +E exc ), where S n is the neutron 
separation energy in the parent nucleus A. The energy in 
parentheses is just the binding energy of the neutron that 
decayed. Then, all other things being equal, extremely 
large values of E exc will be suppressed by the decrease 
of the phase space of the 3 neutrinos. This will tend to 
favor the excitation energies of the states in the daughter 
being similar to the binding energy of the neutron that 
decayed. Of course, this argument also assumes that the 
recoil kinetic energy of the residual nucleus will be neg- 
ligible. 



A. The n — > v v v decay mode in Pb. 

As discussed above, the signature of this nucleon decay 
would be given by the specific properties of the residual 
nuclei resulting from the decay. In this case, 208 Pb would 
first become 207 Pb. What would happen next would de- 
pend on the excitation energy of the states in 207 Pb that 
were populated compared to the one-neutron emission 
threshold, the two-neutron emission threshold, etc., in 
207 Pb. Subsequent neutron emissions will be much more 
rapid than electromagnetic de-excitations, so subsequent 
neutron emissions would occur instead of de-excitations 
of the daughter nuclei as long as the states populated 
were above the neutron-decay threshold. The excita- 
tion energy in the post-decay 207 Pb will depend on the 
binding energy of the nucleon that decayed. This en- 
ergy would be expected to be as much as several tens of 
MeV in lead, but would be weighted toward lower values 



by the preponderance of higher spin nuclear orbits, hence 
higher occupation numbers, near the Fermi surface. Thus 
a reasonable range to assume for the excitation energy in 
207 Pb might be an asymmetric distribution ranging from 
essentially zero (if a valence neutron decayed) to as much 
as 20 MeV, with an even higher-energy tail of the distri- 
bution extending to several more tens of MeV, but with 
the highest energies suppressed. The one-neutron (two- 
neutron) emission threshold in 207 Pb is 6.74 MeV (14.83 
MeV) . Decays of neutrons in the highest-energy occupied 
shells would presumably primarily populate states below 
the one-neutron-emission threshold in 207 Pb. Although 
these states would decay by emitting 7-rays, 207 Pb is sta- 
ble. Thus the 7-rays would be the only signature of the 
neutron decays to the low-lying states, providing a less- 
than-compelling signature of nucleon decay. However, a 
fairly large fraction of the neutron decays in 208 Pb would 
produce at least one neutron emission. These states 
would produce 207_J 'Pb plus j neutrons, with j being at 
least one. 

The branchings into the particle and 7-ray emission 
channels at a given excitation energy of the daughter nu- 
cleus 207 Pb have been calculated. The relevant trans- 
mission coefficients were determined with the same in- 
puts used in the NON-SMOKER statistical model code 
[l7l ITsf . which is often used for astrophysical calcula- 
tions, and has been found to be accurate over a wide 
mass range 0, |2jj . In addition to single-particle emis- 
sion, two-particle emission can be calculated, specifically 
for the case at hand, the two-neutron emission. For two- 
particle emission each transition was followed from a level 
with given spin and parity in 207 Pb to a definite level with 
given spin and parity in 206 Pb, then probabilities for sub- 
sequent neutron emission were determined by summing 
over neutron emission to all possible final states in 205 Pb. 
Up to 20 low-lying experimentally known levels were used 
in each nucleus involved, and a theoretical level density 
[20| was employed above the last known state. The rel- 
ative probability of one- and two-neutron emission as a 
function of excitation energy in 207 Pb is shown in Fig.^ 
Spins from ~ to ^ and both parities were considered in 
207 Pb. The transitions to the states of different spin were 
weighted only with 2J + 1. In order to obtain the cor- 
rect decay probability, the calculated probabilities have 
to be folded with the function that describes how the 
207 Pb levels of a given excitation energy are populated 
in the primary decay event, as described above. The ba- 
sic features of the present results will remain although 
the actual distribution of strength would be expected to 
be distorted somewhat due to the folding with the pop- 
ulation derived from vvv decay. 

As can be seen from Fig. ^ the neutron channel dom- 
inates the possible decays above the neutron emission 
threshold. Only up to about 9 MeV is the decay gov- 
erned by photon emission. Although the proton emission 
threshold is at 7.5 MeV, there is no significant proton 
emission up to 40 MeV in excitation due to the high 
Coulomb barrier; it is still only 20% of all decays at 
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FIG. 1: The relative probability for one- and two-neutron 
emission as a function of of excitation energy in 207 Pb. 



48 MeV of excitation energy. Although the two-neutron 
emission channel opens at 14.83 MeV, it remains insignif- 
icant up to 30 MeV in excitation energy reaching 40% of 
all transitions at 48 MeV. It should be noted that the 
solid curve in Fig. ^ contains the sum of one- and two- 
neutron emission. Therefore, emissions of single and dou- 
ble neutrons contribute equally at the highest calculated 
energy, each comprising 40% of all decay possibilities. 
Thus the population of excitation energies close to the 
neutron separation energy can be probed by measuring 
7-rays, whereas higher excitation energies are accessible 
by two-neutron decay. 

The lowest lying 82 neutrons in 208 Pb would be suf- 
ficiently strongly bound that their decay would be very 
likely to populate states at least above the one-neutron 
emission threshold. Some of the remaining 44 neutrons 
would also be likely to produce daughter states above 
that threshold. Their ordering is h g / 2 , £7/21 ^5/2) P3/2> 
ii3/2i an d Pi/2j an d their approximate binding energies 
are 12, 11, 8, 8, 8, and 7 MeV respectively [2l|. It seems 
reasonable to assume that decay of the neutrons in the 
last three orbitals would populate states below the one- 
neutron emission threshold. Thus 106 of the 126 neutrons 
in 208 Pb would be expected to populate states above the 
one-neutron emission threshold in 207 Pb. An additional 
concern, though, is that the decays of neutrons in the 
most tightly bound states, i.e., those dominated by the 
lsi/2, IP3/2, lpi/2, ld 5 /2, 2si/ 2 , and ld 3 / 2 orbitals, might 
result in states in 207 Pb that would decay by proton or 
a-particle emission. Their binding energies are 40, 35, 
35, 32, 31, and 29 MeV respectively |2lJ. However, our 
emission probability calculations described above suggest 
that only the states having ls!/ 2 neutrons would be ex- 
pected to exhibit such decays; we have subtracted the 
two neutrons in that orbital from consideration. Thus 
104 neutrons in 208 Pb can decay to produce detectable 
signatures. 

The upshot is that this nucleon decay mode would have 
a high probability, roughly 104/126 (so that e n (see eq. 1) 



is 104/126=0.825), for producing at least one neutron in 
coincidence with a fairly high-energy 7-ray, but emission 
of several neutrons in coincidence with a 7-ray would 
also have a relatively high probability. We have assumed 
50% for the sum of the probabilities of events in which at 
least one neutron and a detectable (> 3 MeV) 7-ray are 
emitted. We note that states above 4 MeV of excitation 
in 207 Pb have a high probability, at least 50% 0], of 
producing a 7-ray of at least 3 MeV in energy. 

We assume 1 kT of lead perchlorate admixed with 20% 
water (0.41 kT of lead); that has 5.9xl0 29 208 Pb nuclei, 
or 6.1xl0 31 neutrons in 208 Pb that would be expected to 
populate states above the one-neutron emission threshold 
of 207 Pb and decay subsequently by neutron emission. If 
the fraction of the 104 neutrons assumed that decay by 
single or multiple neutron emission is 100%, the probabil- 
ity for detection of at least one neutron is 50% (it will be 
considerably higher for multiple neutron emission), the 
probability for producing a detectable 7-ray is 50% and 
its detection efficiency is 50%, then the probability that 
a neutron decay event from one of the 104 neutrons as- 
sumed to be detectable will be observed is 12.5%. If the 
lifetime for this decay branch is 10 30 years then, under 
these circumstances, 7-8 events per year would be ob- 
served. Note that we have assumed that only the 208 Pb 
would contribute; it is likely that the other lead isotopes 
would produce similarly detectable decays, producing up 
to a factor of 2 enhancement. Furthermore, the assumed 
efficiencies are conservative. Thus this mode could have 
its lifetime extended by searching for these decays in lead, 
provided the background events could be managed. 

The most obvious background signal that has several 
neutrons that wouldn't be vetoed by a cosmic ray shield 
would involve production of those neutrons by an ener- 
getic cosmic ray in the surrounding rock. One might then 
get a 7-ray if one of the neutrons could inelastically excite 
a nucleus. However, this background could be eliminated 
easily, as the LPC detector could be surrounded by mod- 
erator, so that none of the neutrons getting into it would 
have enough energy to inelastically excite a nucleus. 

A more serious background would result from neutri- 
nos produced in the Earth's upper atmosphere inelasti- 
cally exciting 208 Pb via the neutral-current interaction 
to levels from which it could emit one or more neutrons, 
going to nuclei that might themselves emit more neu- 
trons. This background would be impossible to reject 
on an event-by-event basis. However, its yield could be 
estimated as a function of the energy of the incident neu- 
trino by measuring the number of charged-current inter- 
actions as a function of energy (as determined by the 
energy of the recoiling lepton) that would produce sim- 
ilar numbers of neutrons in coincidence with the 7-ray. 
One could then infer the yield from the neutral-current 
interactions on the basis of the relative magnitudes of the 
two types of cross sections and the energy distribution of 
the charged-current interactions. 
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FIG. 2: The relative probability for one- and two-neutron 
emission and one-proton emission as a function of excitation 
energy in 34 CI. 



B. The n — » v v v decay mode in CI. 

The abundance of 35 C1 is 75% of natural CI, and there 
are 2 CI atoms per Pb atom in LPC. Thus it is also useful 
to see if nucleon decay in 35 CI might produce a definitive 
signature. Neutron decay in 35 CI would sometimes pro- 
duce 34 C1 in a highly excited state. However, the one- 
proton emission threshold for 34 C1 is at 5.14 MeV, far 
below the one-neutron emission threshold at 11.51 MeV. 
Indeed, there are apparently no bound levels in 3 CI even 
close to 11.51 MeV. Thus the result of neutron decay to 
34 C1 will be 7-ray decays to the ground state of 34 C1. 
Note, though, that 34 C1 has an isomeric state at 0.146 
MeV; so roughly half of the energetic 7-rays would go to 
that state. 

The relative branchings into different decay channels 
of 34 C1 were calculated using the same methods as de- 
scribed above for 207 Pb. The results are shown in Fig. [21 
The situation for 34 C1 is different than for 207 Pb, how- 
ever, because of the lower Coulomb barrier and a proton 
separation energy that is lower than the neutron separa- 
tion energy. Neutron emission is relatively unimportant 
at all calculated excitation energies. The excited nucleus 
34 C1 will de-excite via 7-transitions for excitation ener- 
gies up to about 9.5 MeV, where proton emission ac- 
counts for 50% of all emissions. Proton emission quickly 
rises to nearly 100% by 11.9 MeV and the falls off at high 
energies due to increased a emission. Two-neutron emis- 
sion is completely negligible at all calculated energies. 
As with Pb, a quantitative description of the emission 
requires the knowledge of the population of the excited 
states by the decay, but would not be expected to differ 
qualitatively from the results shown here. 

The ground state of 34 C1, however, /3-decays with a 
half-life of 1.53 s. Thus a coincidence between the en- 
ergetic 7-rays (assumed detection efficiency = 50%) and 
the /3 + -decay would identify a candidate event for neu- 
tron decay in 35 CI. The end point energy of the ft + is 4.47 



MeV, so most of the ft + s will be produced with sufficient 
energy to be detected from their Cherenkov radiation (as- 
sumed detection efficiency = 70%). The isomeric state 
has a much longer half- life (32 m), so we have assumed 
that would not produce a useful coincidence for deter- 
mining nucleon decay (reducing the useful event yield by 
half). None the less, nucleon decay in 35 C1 could be iden- 
tified by observing ft and 7s in coincidence, with enough 
events detected to confirm the half- life of 34 CI. The rela- 
tively long half-life of 34 C1 would demand that the mate- 
rials used in OMNIS LPC modules be as pure as possible 
to minimize accidental coincidences. 



Again assuming 1 kT of lead perchlorate admixed with 
20% water, there will be 1.8xl0 30 35 C1 nuclei. Of the 18 
neutrons in 35 C1 those in the lsi/2, IP3/2, and lpi/2 or- 
bits (binding energies = 36, 27, and 23 MeV respectively 
|2l|'). would probably decay to sufficiently highly-excited 
states in 34 CI that they would emit a proton, neutron, 
or a-particle, and not end up in 34 C1. The remaining 
three orbitals, S1/2, and d^/o have binding energies 
of roughly 17, 14, and 10 MeV 21]. The one-neutron 
separation energy in 35 C1 is 12.64 MeV; removal of a va- 
lence neutron, presumably from a d3/2 dominated state, 
would require that much energy and would tend to pop- 
ulate the ground state of 34 C1. Thus, presumably, decay 
of neutrons in 35 CI states having strong ld 5 / 2 or 2si/2 
neutron configurations would tend to populate states in 
34 C1 of sufficient excitation energy to emit a detectable, 
> 3 MeV, 7-ray. Thus, about 8 neutrons per 35 CI nu- 
cleus, or 1.4xl0 31 neutrons in the 35 C1, (the product N n 
e„ in eq. 1) would have appropriate binding energy to 
decay to states in 34 CI that could produce a 7-ray above 
3 MeV in energy, so would have a high probability of 
being detected. If the detection efficiency for the ft is as- 
sumed to be 70%, and the lifetime for this decay process 
is 10 30 years, one would expect to observe 3-4 events per 
year from neutron decay in 35 CI. Although this decay 
mode would not produce as strong a limit on nucleon 
decay as would the lead, its signature would be consider- 
ably more definitive than that from the lead, because its 
backgrounds are so much more readily identifiable. 



The primary background for decay of a neutron in 
35 CI comes from neutrinos produced in the Earth's at- 
mosphere, but such events would not be able to simulate 
nucleon decay events. Although once the neutron was 
knocked out of the 35 CI, the resulting 7 — ft coincidence 
would be the same as for the neutron decay event, the 
atmospheric neutrino would emit a neutron, which would 
be detected to give its characteristic 8.6 MeV 7-ray. This 
is a considerably greater energy than could be produced 
by any of the 7-rays resulting from the de-excitations in 
34 C1; thus that neutron could be used to veto this type 
of background event. 
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C. Monte-Carlo simulations of detectors to 
determine detection efficiencies 

The limit that can be achieved from any process such 
as the decays described above clearly depends on the de- 
tection efficiency. Thus we have run simulations of the 
events from decays of 35 CI based on the GEANT de- 
tector simulation software. As planned, the LPC detector 
will consist of multiple independent cylindrical modules 
having radii of 2 m and heights, which are adjustable in 
the simulations, ranging from 1 to 2 m. The GEANT 
simulations were based on a single module. The inte- 
rior of the cylinder was assumed to be viewed by photo- 
multiplier tubes (PMTs) at both ends. Current design 
plans are to use an LPC solution of concentration, also 
adjustable in the simulations, of 75% to 80% LPC by 
weight. 

From the decay of a neutron in 35 CI, GEANT gener- 
ated both the 7-rays from the de-excitation of the 34 CI 
nucleus and the resultant (3 + from the decay of the 34 C1. 
The initial starting state for 34 CI was randomly chosen 
from the energy levels that were below the proton sep- 
aration energy. 7-rays were then generated to simulate 
the de-excitation of the 34 C1 nucleus from the selected 
starting state down to either the ground state or the iso- 
mer. The number and energies of the generated 7-rays 
were based on data on the levels and 7-rays for the 34 C1 
nucleus. Events that decayed to the isomeric state were 
rejected. The fraction of the events that de-excite down 
to the ground state was found to be 57.6 ± 0.2%. 

The 34 CI decay time is chosen by randomly selecting 
whether a decay will occur during the 1 ms time unit, 
and was repeated until a decay occurred. The chance 
for a decay within a time unit was based on the half- 
life of 34 CI. If the decay time took longer than 40 s, the 
total accumulated time was set to zero and the decay test 
continued. 

The kinematics for the (3 + were generated after the 7- 
rays were generated and tracked. The initial energy of 
the (3 + was selected from a lookup table and the decay 
time was added to the initial time of flight. A Gaussian 
profile was used for the initial energy distribution of the 
positron. The parameters for the Gaussian were derived 
from a fit to the /3 + energy spectrum calculated using the 
RADLST program ||. 

Upon entering a PMT, a photon's energy and time 
of flight were stored. As f3 + s and /3~s have a large 
RMS multiple scattering in LPC, any patterns in the 
Cherenkov radiation are destroyed. However, an event 
can be identified by the number of PMTs that have fired 
in a localized area within a specified time window. The 
requirement imposed to identify 7-rays and /3 + s required 
three PMTs to fire in a localized cluster within a 20 ns 
window. 

To estimate the efficiency for detecting neutron decay 



in the 35 C1 in LPC, 10,000 events were generated in a 100 
cm deep tank with an 80% solution of LPC by weight. 
The attenuation length of the LPC was taken to be 4.2 
m The mean efficiency for detecting the 7-rays from 
the de-excitation of the 34 C1 nucleus over the volume of 
the detector was found to be 57.0 ± 0.8%. For detecting 
the /3 + s, the mean efficiency was found to be 75.2 ±0.9%, 
making the mean efficiency for detecting a coincidence 
between the 7s and the /3s to be 42.6±0.8%. For a 200 cm 
deep tank with an 80% LPC solution, the mean detection 
efficiency for 7s, /3 + s, and coincidences was found to be 
42.6 ± 0.7%, 60.6 ± 0.9%, and 25.6 ± 0.5% respectively. 
Thus the efficiencies assumed above are consistent with 
those determined from the Monte-Carlo simulations. For 
the 200 cm deep tank, the neutron detection efficiency 
was found to be 87.6 ± 0.9%. 

The position of a nucleon decay event can be deduced 
from the time difference At between hits at the two ends 
of the tank. At can be defined in a variety of ways, e.g. 
by taking it to be the time difference between the peaks of 
the signals, either from the 7s or the /3 + s, from the two 
sides of the module, or by attempting to average early 
hits to utilize the leading edges of the signals on the two 
sides to do the timing. The first approach can be con- 
fused somewhat by the arrival of photons scattered from 
the opposite side, whereas the latter can be complicated 
by limited statistics. Either approach seems to allow lo- 
calization of events to about ±10 cm. Spatial localization 
is important, as it allows use of a position dependent neu- 
tron detection efficiency, so can allow greater emphasis 
on events that occur in the center of the detector, which 
will have a greater efficiency for detecting the veto neu- 
tron, than those at the edges, which will have a lesser 
efficiency. 



III. CONCLUSIONS 

Given that the suggested searches for nucleon decay 
would be conducted in a supernova neutrino detector, 
these searches would have a long time, probably more 
than 20 years, to run. Thus the suggested experiments 
have the potential to improve greatly the existing limits 
on the n — > v + v + v decay mode. 
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